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Abstract 

We present a strategy for a lattice evaluation of the neutron electric dipole moment 
induced by the strong CP violating term of the QCD Lagrangian. Our strategy is 
based on the standard definition of the electric dipole moment, involving the charge 
density operator Jo, in case of three flavors with non-degenerate masses. We present 
a complete diagrammatic analysis showing how the axial chiral Ward identities can 
be used to replace the topological charge operator with the flavor-singlet pseudoscalar 
density Pg. Our final result is characterized only by disconnected diagrams, where the 
disconnected part can be either the single insertion of Pg or the separate insertions of 
both Pg and Jo- The applicability of our strategy to the case of lattice formulations 
that explicitly break chiral symmetry, like the Wilson and Clover actions, is discussed. 
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1 Introduction 



A non-vanishing value of the electric dipole moment (EDM) of the neutron is a direct 
manifestation of the breakdown of both parity and time reversal symmetries. Within the 
Standard Model and its possible extensions parity and time reversal can be violated both in 
the electroweak and strong sectors. In this study we focus on the neutron EDM induced by 
the so-called #-term of the QCD Lagrangian [T|, given in the Euclidean space by 

g 2 ~ g 2 

where G° is the gluon field strength, c the color octet index (c = 1, ...,8) and 9 a dimen- 
sionless parameter. 

The present experimental upper limit on the neutron EDM, = \d^\ < 6.3 • 10~ 26 (e ■ 
cm) at 90% confidence level [2|, which corresponds to a severe bound on the magnitude 
of 9. Indeed, using the available theoretical estimates from Refs. Ej, one has « 
3 • 10~ 16 \9\ (e • cm) leading to \9\ < 2 • 10~ 10 . The smallness of the parameter 9 is usually 
referred to as the strong CP problem. 

Available estimates of the relevant matrix element however are based on phenomenolog- 
ical models, as the MIT bag model of Ref. [3j or as the effective nN chiral Lagrangian of 
Ref. j3j. Estimates relying on non-perturbative methods based on the fundamental theory, 
like lattice QCD, are still missing to date. A first attempt to calculate d^ with lattice QCD 
was done in Ref. [S] long time ago. There, using the axial anomaly (as suggested in Ref. 
the ^-term (JIJ was replaced by the flavor-singlet pseudoscalar density Pg, obtaining 

% % r — " 

C'e = Tj-9 mP s -> --9 m u^u + d^d + sj 5 s , (2) 
/ 3 J 

where YrT 1 = (m' 1 + m^ 1 + m~ 1 )/3 in case of three flavors with non-degenerate masses. 
Then, d^ was evaluated by extracting the spin- up and spin-down neutron masses from the 
two-point functions obtained by adding to the action both C'e and a term corresponding to 
a constant-background electric field oriented in a fixed spatial direction. Since 9 is expected 
to be quite small, it is enough to consider a single insertion of C'e in the two-point functions. 
The latter can be divided into connected and disconnected insertions, when the operator (J2J) 
is contracted in a valence-quark line and in a virtual quark loop, respectively, as pictorially 
represented in Fig. [TJ 

In Ref. jS] a lattice calculation using Wilson fermions and including only connected 
diagrams was attempted and a no n- vanishing result for d^ was found. In Ref. jH] it was 
shown however that the connected diagrams cannot contribute to the neutron EDM and 
therefore the non-zero result of Ref. [3] arises from lattice artifacts of 0(a). Thus a lattice 
QCD estimate of the neutron EDM is not yet available. 

It should be pointed out that the approach of Refs. E] is based on the technique of 
introducing an external constant-background field. This approach has been abandoned since 
modern lattice QCD simulations allow to evaluate directly the relevant matrix elements of 
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Figure 1: Connected (a) and disconnected (b) insertions of the flavor-singlet pseudoscalar 
density operator. The insertions are represented by open dots, while solid lines are quark 
propagators. The hatched ovals are the operators which create and destroy the neutron. 
Gluon lines as well as extra quark loops are not shown. 

the electromagnetic (e.m.) current operator. Moreover, in the proof of Ref. [6J it is assumed 
that the connected insertions of the singlet and non-singlet pseudoscalar densities are the 
same. Since this is not the case as we have explicitly checked, a more careful diagrammatic 
analysis of the insertions of the pseudoscalar operators appears to be worthwhile. 

In this Letter we present an alternative way to compute the neutron EDM on the 
lattice, which avoids the introduction of the background electric field thanks to the use of 
the standard definition of the neutron EDM 



where Jo is the time component of the e.m. current operator. If the #-term is neglected 



where \N) is a shorthand for \N) e=0 . 

As well known [7] , the direct calculation of Eq. (J3J) is a difficult task due to the presence 
of the topological charge operator. Therefore it is of interest to investigate the possibility 
of replacing the insertion of the topological charge with the disconnected insertions of the 
flavor- singlet pseudoscalar density in the presence of the charge density operator J . Since 
the insertions of the latter can generate disconnected diagrams, our analysis represents also 
an extension of the previous one made in Ref. (H]. 

The plan of the paper is as follows. In Section El we briefly recall the non-singlet and 
singlet continuum axial Ward Identities (WPs) which are relevant in this work. In Section El 
we present a complete diagrammatic analysis, showing that the topological charge can be 
replaced by the flavor- singlet pseudoscalar density in Eq. (0J. Our final result is characterized 
only by disconnected diagrams, where the disconnected part can be either the single insertion 
of Ps or the separate insertions of both Ps and Jo. The applicability of our strategy to the 
case of lattice formulations that explicitly break chiral symmetry, like the Wilson and Clover 
actions, is discussed. Our conclusions are summarized in Section |U 




(3) 



in the QCD Lagrangian, then the moment d N identically vanishes. Treating Cg as a pertur- 
bation at first order, one has 
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2 Axial Ward identities 



Axial WPs can be obtained in the usual way starting from the vacuum expectation value of 
a generic operator O. In Euclidean space one has 



(0(x h ...,x n )) = — I d[G\d[^\d[^\ 0(x 1 ,...,x n ) e 



(5) 



where S is the QCD action without the #-term and Z = J d[G]d[ip]d[tp] e~ s . From now on 
the notation (...) refers to vacuum expectation values evaluated without the #-term (JTJ in 
the action. Performing the non-singlet axial rotations of the quark fields ip(x) and ip(x) 



^(x)^[l+ia a (x)X a l5 ]il;(x) , 



where A a are the usual SU (3) flavor matrices (a 



(O 



5S 



) = (- 



50 



5(ia a (x)) 5(ia a (x)) 



8), one gets 



(6) 



(7) 



with 



(O 



5S 



X a 



5(ia a (x)) 



) = -d,(OA a (x)) + (O ^(x) M, — 75 ^(x) ) 



where M is the SU(3) mass matrix and A a (x) the non-singlet axial current. From now on 
all the relevant operators and fields should be considered as renormalized quantities in the 
continuum. The case of the lattice regularization will be addressed at the end of Section El 
Performing the flavor-singlet rotations 



one gets 



ip(x) 
4){x) 

(O 



1 + ia (x)'~f5 ip(x) 
i>(x) 1 + za°(a;)75 



ss 



:) = <■ 



50 



5(ia°(x))' "5(ia°(x)) 



> 



(9) 



(10) 



where 



(O 



5S 



6{ia o {x)) ) = -WOA%x)) + (O ^(x) {M, A }75 1>(x) ) + 2N f ^- 2 (OGG(x)} , (11) 



with Ao being the SU(3) identity matrix and A^(x) the (renormalized) flavor-singlet axial 
current. In Eq. (|TT|) the anomalous term proportional to the (renormalized) GG operator 
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has been taken into account and g is the renormalized strong coupling. Note that in the non- 
singlet channel all the terms appearing in the r.h.s of Eqs. (jHJ) are separately renormalization- 
group invariant, while in the singlet case [Eq. a divergent mixing between d^A^ and 

g 2 GG appears. 

Following the trick suggested in Ref. |3j, Eqs. (fTUJ) and (JHJ) can be used to replace GG(x) 
in Eq. (jlj) with the pseudoscalar density Ps{x). By integrating the WI over the whole space- 
time, the divergence of the axial current d^A^, appearing in the r.h.s. of Eq. vanishes 
identically a , whereas the flavor- singlet contact terms given by the variation of O in the r.h.s. 
of Eq. (|T0|) . may give a non-null contribution. 

Without loss of generality we can limit ourselves to the simple S*?7(2)-symmetric case, 
where m u = m^ = m, and m s ^ m. Since 



{M, A } = 2m 



1 + 



An 



Vs 

V3 



(12) 



where r] s = (m s — m)/m, the singlet WI involves both the flavor-singlet pseudoscalar 
density Ps and the non-singlet operator P 8 . To eliminate the latter we have to combine the 
singlet WI (fTT]) and the non-singlet WI (JHJ) with a = 8. By exploiting the relation 



M, 



A* 



2m 



1 + 



277, 



A« 
2 



Vs 



3V3 



Ar 



(13) 



one finds 



{M, A } + ^ 1 —— Im, - 



2m\r 



with m 



m 



(14) 



~ x + m 1 + m s 1 )/3. Equation (fTlj) implies 



2N f 



32tt 2 



(i 4 x (OGG(x)) = (5[0)) - 2m / d 4 x (OP s (x)) 



with 



(S[0\) = {S s [0}) + 



2Vs 



1 



1 + 2r) s /3 



(15) 



(16) 



where 5s [O] and 5%[0] stand for the variations of O integrated over the whole space-time and 
corresponding to the singlet and octet chiral rotations [see Eqs. (JHJ and ©], respectively. 
In the next Section we present a complete diagrammatic analysis of the r.h.s. of Eq. (|15|) . 

a This is true in absence of massless Goldstone-boson poles, which certainly occurs away from the chiral 
limit. In case of the singlet channel the space-time integration of d^A^ vanishes also in the chiral limit of 
full QCD, when the 7/-meson remains massive because of the gluon anomaly. Let us remind however that 
in the chiral limit the 9-term Q can be eliminated by a chiral rotation of the quark fields and therefore in 
that limit <ijv must vanish. 
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3 Neutron EDM and the disconnected insertions of the 



pseudoscalar density 



In the flavor- singlet channel the operator O of interest for Eq. (jlj is given by 



O = N a (z) J (y) N p (0) 



(17) 



where Jo(y) is the e.m. charge density 



My) = e « u (yho u (y) + e d d{y)^ Q d{y) + e s s 



(yhos(y) 



(18) 



and N a (p) are the interpolating fields of the neutron jHj 



N a = -U^ (d T ^ 5 u k ) 



(19) 



with k being color indexes, a and (3 Dirac spinor indexes, and C the charge conjugation 
operator. As well known, standard procedures in lattice QCD calculations allow to extract 
the relevant matrix element using the interpolating fields (|19j) at enough large values of the 
source (sink) time distances. 

Let us now consider the chiral variations 6s[0) and Ss[0] of the operator (|T7j). which can 
be concisely written in the form 



5j[o\ = sANa&z) My) n p (o) + N a { z ) 5j[M(y) mo) + N a {z) My) (o) (20) 



where j stands for j = S or j = 8. Through the package FORM jS] we have carried out 
all the Wick contractions involved in the variations (|2(Jj) as well as in the quantity OPs(x). 
The Feynman diagrams corresponding to the chiral variations 6s[0] and S$ [O] cannot be 
directly compared with those corresponding to OPs(x), because the former contain one 
quark propagator less than the latter. In order to make the comparison possible we make 
use of the non-singlet axial WPs (J7J) and (jHJ) with the operator O given by O = ij){w) i(j(w'). 
After integration over the whole space-time one obtains 



where Si is the propagator of the i-th quark. The systematic use of Eq. (|2T|) in the Feynman 
diagrams generated by the chiral variations of the operator O allows to gain one more 
propagator and one more integration in a such a way that the diagrams arising from the 
chiral variations can be now directly compared with the Feynman diagrams generated by 
the operator OPs(x). The latter can be divided into connected and disconnected diagrams, 
whose possible typologies are sketched in Fig. |21 The chiral variations cancel out all the 
connected diagrams [see Fig. I2te)], and therefore we confirm the final result of Ref. that 
connected diagrams do not contribute to the neutron EDM. However, the chiral variations 
contain also disconnected diagrams of the type (c) and (d) which cancel out the corresponding 



Si(w, w') 75 + 75 Si(w, w') = 2mi / d 4 x Si(w, x) 75 Si(x, w') 



(21) 
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(a) (b) 



(c) (d) (e) 

Figure 2: Possible typologies of disconnected diagrams [(a) - (d)] and an example of a con- 
nected diagram [(e)], involved in OPs(x). The meaning of the vertices and lines is the same 
as in Fig. QJ Open dots and full squares represent the insertion 0/75 (Ps) and 70 (Jo), 
respectively. 

diagrams present in OPs(x). While diagrams of the type (c) vanish identically, as it can 
be checked using Eq. (}2"T]) and the relation {75,70} = 0, the diagrams of the type (d) are 
non-vanishing in case of 5'[/(3)-symmetry breaking. Therefore, our general result is 

2 

2iV/ db / * X (®GG(x)) = -2m J d*x { [(OP s (x))) disc _ {a) + i(OP s (x))] disc . (b) } , (22) 

which allows to replace in Eq. (0J) the topological charge with disconnected insertions of the 
pseudoscalar density. A pictorial representation of all the explicit diagrams contributing to 
[(OPs(x))} dlsc {a) and [(OP s (x))} disc {b) is shown in Fig. 

Note that: i) in the SU(3) symmetric limit the diagrams of type (d) vanish, because the 
disconnected insertions of the charge density Jo are zero due to the relation e u + e^, + e s = 0; 
ii) the r.h.s. of Eq. (J2~2"j) is proportional to m. Since fn = when at least one of the quark 
masses is zero, the insertion of the topological charge has no effect in such a limit and the 
neutron EDM is vanishing; iii) the only insertions of the singlet pseudoscalar density, 
which are left in the final result ()22|). are those disconnected diagrams which are related to 
the operator GG via the anomaly. This is why the diagrams of the types (c) and (d) of Fig. El 
do not contribute to the neutron EDM. 

The final result of Ref. [H] can be recovered from Eq. (|2*2"|) simply by considering the SU(3) 
symmetric limit and by dropping the charge density Jo in the operator O [see Eq. (j!7j) ]. In 
this case [{OPs(x))] disc , b -, = and the disconnected contributions [(OPs(x))] disc , a -> are those 
illustrated in Fig. HI 

Before closing this Section, we show the applicability of our results to the case of lattice 
formulations, that break explicitly chiral symmetry, as the Wilson or Clover actions. To do 
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K 0P s)] disc . (a) =- N c 



(N f -1) 



A 



+ 



i e u+ e d) 



+ e d 



+ e d 



-(H)- •*-(&$■ 



I(° p s)] disc , b) = +N - 



(N f -l) 



+ 







+ e.l ) 



Figure 3: Disconnected diagrams contributing to the r.h.s. of Eq. \2fy) . The notation is the 
same as in Fig. HJ but solid lines are now u- and d-quark propagators, while dashed lines are 
s-quark propagators. 



that we consider the works of Refs. ^01 ITT] , whose central result is that even within the 
above-mentioned formulations the continuum limit can be reached in such a way that the 
continuum WPs are respected. 

When the functional integral is regularized on the lattice using Wilson fermions the 
effect of a non-singlet axial rotation leads to the following identity: 



A,(OA a (x)} + (O fax) I M , 75 fa)) + (OX a (x)) = ( 



S[0] 



S(ia a (x)) 



) 



(23) 



where and ip are a lattice version of the (non-singlet) bare axial current and quark field, 
respectively, and Mq is the bare SU(3) quark mass matrix. In Eq. (|2*3j) A M is the backward 
derivative on the lattice, while the operator X a (x) is the chiral variation of the Wilson term. 
The latter is known to mix with operators of the same and lower dimensionality, namely 



(OX a (x)) 



1 - z 



(NS) 



A,<04°(*)> - (O fax) |M (ArS) ,y| 7 5 fax)) - (OX a (x)) (24) 
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Figure 4: Disconnected insertions of the flavor- singlet pseudoscalar density when the SU(3) 
symmetric limit is considered and the operator O [Eq. fli7| )/ does not contain the charge 
density operator J . The meaning of the vertices and lines is the same as in Fig. [3J 



where Z 



vs ^ and M^ NS ^ are mixing coefficients, and X a (x) is an operator whose matrix 



elements vanish in the continuum limit. Substituting Eq. 



-A M <0 



r(NS) 



K(x) > + (o 



into Eq. (|2"3j) one gets 
8[0\ 



) = { 



S(ia a (x)) 



) + 0(a)(25) 



where M^^ = [M 



M iNS) yz 



(TVS) 



with Z^p S ^ being the renormalization constant of 



the non-singlet pseudoscalar density. Therefore, one can identify the continuum limit of 
Z^ S ' A^ix) and Z^p fS ^tfj(x) ^M^ NS \ A a /2 j 75 i/j(x) with the renormalized axial current A a ^(x) 
and ijj(x) {M^ ns \ A a /2| 75 VK^)) respectively, appearing in Eq. (|SJ). In this way the non- 
singlet continuum WI is recovered with a renormalized quark mass matrix given byM^). 
In the singlet case the chiral variation of the Wilson term, X s , obeys the relation 



(OX s (x)) 



,(S) 



A^(OA^(x)) + counterterms — (O ■?/>( 



+ 



2N, 
32tt s 



;Z {qIgg) 



x 



g (OGG(x)} — counterterms — (OX (x)) 



(26) 



A and Z^2 G q^ are mixing coefficients, GG(x) is a lattice version of the corresponding 



where Z 

classical operator and the counterterms make separately finite the operators A^A^x) and 

#0 GG(x). A complete one-loop analysis of the renormalization constants appearing in 
Eqs. (I24|) and (J26j) for both the Wilson and Clover actions will be reported elsewhere 
The main difference between the singlet and non-singlet channels is that beyond one loop the 

divergence of the axial current Z^f* A^A^x) and the anomaly Z^ 2Q ^ g\ GG(x) are not yet 
finite operators, and a logarithmically divergent mixing, represented by the counterterms in 



Eq. (j2T?j) . is required in order to end up with finite operators (see, e.g., Ref. 



Denoting 



such operators by A^A^'^'^x) and g 2 GG 



(ren.) 



respectively, one gets 



- A,(0 [Af> ren -\x)}) + (O \z^$(x) {M^ s \ Ao} 75 Hx) 



+ 2N f ^- 2 (OGG (x)) = ( 



6[0] 



5(ia°(x)) 



) + 0(a) 



(27) 



where = [Mo — M }/Zp with Zp being the renormalization constant of the singlet 
pseudoscalar density. Therefore, in analogy with the non-singlet case one can identify the 

/m - r i * ' — ~( ren -) 

continuum limit of the operators Jjf ,ren ^ (x) , Z P ip(x) lM ( - s \\ >'j 5 tjj(x) and GG (x) 

with the renormalized axial current A^(x), tp(x) |m^ 5 \ Aoj 75 ip(x) and gluon anomaly 

GG(x), respectively. In this way the singlet continuum WI is recovered with a mass 

M^. Although one has Z { p ] ^ Z P NS) and M (S) ^ M (NS \ it has been shown 14 that 

[M - M (5) ]/^p } = l M o ~ M^^/Z^, which means = M^ NS "> = M. In this way 

the continuum WPs are preserved. In case of S*?7(3)-symmetry breaking the equivalence 
between M (5) and M (JVS) is essential for the applicability of our strategy leading to Eq. (|22j). 

4 Conclusions 

In this paper we have presented a strategy for evaluating the neutron EDM on the lattice 
induced by the strong CP violating term of the QCD Lagrangian. Our strategy is based 
on the definition of the neutron EDM given by Eq. (0J), which involves the insertion of the 
topological charge in the presence of the charge density operator Jo. In case of three flavors 
with non-degenerate masses we have presented a complete diagrammatic analysis showing 
how the axial chiral Ward identities can be used to replace the topological charge operator 
with the flavor-singlet pseudoscalar density P5. Our final result ()22|) is characterized only 
by disconnected diagrams, where the disconnected part can be either the single insertion of 
Ps or the separate insertions of both P$ and Jo- In this way we confirm the final outcome 
of Ref. that connected insertions of the pseudoscalar density do not contribute to the 
neutron EDM. Finally we have discussed the applicability of our strategy to the case of 
lattice formulations that break explicitly chiral symmetry, like the Wilson and Clover actions. 
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